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I  .  INTRODUCTION 


In  recent  years  rail  guns  have  been  used  to  accelerate  macroscopic  pro- 

1  2 

jectiles  to  velocities  as  high  as  several  km/s.  9  Some  interesting  results 
obtained  in  a  recent  experiment,  to  be  discussed  presently,  have  caused  specu¬ 
lation  that  atmospheric  air  in  the  gun  tube  may  significantly  limit  projectile 
velocity.  In  this  communication  we  wish  to  determine  under  what  conditions 
such  an  effect  might  be  important. 


II.  CALCULATIONS 

The  principle  on  which  the  rail  gun  operates  is  shown  schematically  in 
Fig.  1.  Current  is  conducted  along  the  rails  in  the  direction  indicated  and 
from  the  lower  to  upper  rail  through  an  armature  which  may  be  either  solid  or 
gaseous.  The  resulting  magnetic  field  then  interacts  with  the  armature  cur¬ 
rent,  giving  rise  to  a  force  which  accelerates  the  armature  and  projectile 
(shaded  area)  down  the  gun  tube.  Unless  the  region  of  the  gun  tube  in  front 
of  the  projectile  is  evacuated  prior  to  acceleration,  the  projectile  motion 
will  produce  a  shock  wave  that  exerts  a  retarding  pressure  Pg  on  the  acceler¬ 
ating  system.  We  will  now  solve  the  equations  of  motion  of  the  armature  and 
projectile  to  determine  the  effect  of  Pg.  It  is  assumed  that  the  current  varies 
exponentially  according  to  the  relation 

i  =  iQ  e  t/t0,  (1) 


since  the  current  in  our  rail  gun,  which  has  inductive  energy  storage,  has  been 
successfully  fit  to  such  a  profile.  Results  should  not  be  affected  qualita¬ 
tively  for  other  types  of  variation.  In  the  case  of  a  plasma  armature  it  is 
also  assumed  that  the  entire  system  (armature  and  projectile)  is  accelerated 
at  the  same  rate. 


Under  these  conditions  the  appropriate  equation  of  motion  becomes 


m 


dv 

dt 


.2 
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e-2t/t0 


A 


(2) 


where  m  is  the  total  mass  of  the  armature  and  projectile,  L^  the  inductance 

per  unit  length  of  the  rails,  and  A  the  cross-sectional  area  of  the  gun  tube. 

Now,  for  typical  rail  guns  and  acceleration  times  of  interest,  the  first  term 

on  the  right-hand  side  of  Eq.  (2)  is  large  compared  to  the  second  except  when 

P  is  much  greater  than  atmospheric  pressure.  Consequently,  we  can  approxi- 

^  3 

mate  P  by  a  strong-shock  approximation,  namely, 

^ S.C .  Rashleigh  and  R.A.  Marshall,  J.  Appl.  Phys.  49_,  2540  (1978). 

2 

R.S.  Hawke,  R.L.  Brooks,  F.J.  Deadriek,  J.K.  Saudder,  C.M.  Fowler,  R.S.  Caird, 
and  D.R.  Peterson,  IEEE  Trans.  Magn.  MAG-18_,  821  (1982). 

3 

Y.B.  Zel'dovich  and  Y.P.  Raiser,  Physios  of_  Shock  Waves  and  Hi£h-Tem£ercrture^ 
Hildoodynomic  Phenomena,  Academic,  New  York,  1966}  Chap.  I. 
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ARMATURE 


Figure  1,  Schematic  Drawing  of  Rail  Gun. 


p  v2 

S  -  2  p0 


(3) 


where  y  is  the  ratio  of  the  specific  heat  at  constant  pressure  to  that  at  con- 

-z 

stant  volume  and  where  pQ  is  atmospheric  density,  i.e.,  1.29  kg/m  . 


Substituting  Eq.  (3)  into  Eq .  (2)  we  obtain  a  nonlinear  differential  equa¬ 
tion  for  the  velocity  v.  The  equation  can  be  linearized,  however,  by  the  trans¬ 
formation 


v  = 


2m  du 


(1+Y)P0  A  u  dt  ' 


and  further  simplified  by  letting 


-t/tn 
z  =  zQ  e  0 


where 


zo  = 


r T  •  2  ,,  .  .  2, 1/2 

[L  iQ  (1+y)  p0  A  t0] 


2m 


(4) 


(5) 


(6) 


Physically,  is  a  measure  of  the  importance  of  the  air  in  gun  tube,  its  de¬ 
celerating  effect  increasing  with  increasing  z^ .  We  obtain  after  some  algebra 


d2u  du 


z  u  =  0 


(7) 


a  modified  Bessel  equation.  The  solution  of  Eq .  (7)  consists  of  a  linear  com¬ 
bination  of  the  functions^  Iq(z)  and  Kq(z) .  Using  this  result  for  u  in  Eq. 

(4) ,  applying  the  initial  condition  v(t=0)  =  0,  and  using  formulas^  for  the 
derivatives  of  the  Bessel  functions,  we  find  for  the  time -dependent  velocity 


v  =  vQ(z) 


2z 


lz0  -  * 


Kj(z)  I1(zQ)  -  KjUq)  Ij(z) 
KqU)  IjCZq)  *  KjCZjj)  I0(z) 


Here  v^  is  the  solution  of  the  problem  in  the  absence  of  air  (p 


0 


(8) 

0)  ,  namely. 


j  *  .2 

V0  =-'°  ° 
4m 


2  2 

(1  -  Z  /Zq) 


(9) 


To  determine  when  the  effect  of  the  air  first  becomes  important,  we  can 
(arbitrarily)  set  v/v^  =  0.95  and  solve  Eq.  (8)  for  z  as  a  function  of  z^ . 


4 

B.M.  Murphy,  Ordinary  Pf/ZerentiaZ.  Eguations  and  Their  Solutions^,  van  Nostrand , 
New  York ,  I960,  Chap.  A1 . 

** Handbook  of  Mathematical  Functions ,  edited  by  M.  Abramowitz  and  I.  Stegun, 
National  Bureau  of  Standards,  Washington,  DC,  1964,  Chap.  9. 
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From  the  resulting  value  of  z,  denoted  by  z  ,  the  "characteristic  time"  fol- 
lows  from  the  relation 


*c  "  l0  ln  ' 

For  acceleration  times  less  than  this  value,  the  effects  of  the  air  are  negli¬ 
gible,  while  for  larger  times  they  are  not.  In  Fig.  2  is  plotted  vs. 

from  a  numerical  solution  of  Eq .  (8)  with  v/vq  =  0*95.  For  values  of  some¬ 
what  above  unity,  the  result  is  linear  and  given  by  the  relation 


z  =  zn  -  /07l  5 
c  0 


(ID 


The  validity  of  Eq .  (11)  can  be  demonstrated  from  asymptotic  expansions  of 

the  Bessel  functions.  Since  z  fz  increases  monotonically  with  z  ,  the  charac- 

c  0  u 

teristic  time  t  decreases  with  zA.  Thus,  the  presence  of  the  air  is  likely 
c  u 

to  be  important  for  large  values  of  z^  or  under  acceleration  conditions  corres¬ 
ponding  to  high  currents,  small  masses,  and  long  time  constants. 


We  can  also  determine  from  Eq .  (8)  the  time  at  which  the  maximum  velocity 
occurs.  This  value  is  obtained  when  the  term  in  square  brackets  in  Eq.  (8)  is 
unity  as  can  be  shown  by  setting  the  derivative  equal  to  zero  in  Eq .  (12)  and 
comparing  the  resulting  expression  for  v  with  that  in  Eq .  (8).  Again, we  solve 
numerically  for  z  as  a  function  of  z^  and  denote  the  values  by  .  Results 

are  plotted  in  Fig.  2.  For  zQ  sufficiently  large,  z^  varies  nearly  linearly 

with  Zq  and  satisfies  the  transcendental  relation 


z 

m 


in  (4  z  ) 
v  nr 

2 


(12) 


as  can  be  shown  from  the  appropriate  asymptotic  expansions.  The  time  at  which 
the  maximum  velocity  occurs,  of  course,  is  obtained  from  z^  via  the  analog  to 
Eq.  (10). 

To  demonstrate  further  the  results  of  Eq .  (8)  we  have  plotted  v  and  v^  as 

a  function  of  t  for  data  corresponding  roughly  to  two  recent  experiments.  Eac 

was  carried  out  with  the  same  rail  gun  for  which  appropriate  parameters  were 

A  =  1.61  cm2,  i  =  150  kA,  and  tQ  =  752  ps ;  the  inductance  per  unit  length  L' 

was  assumed  to  be  0.42  pH/m  and  the  parameter  y  was  given  by  1.4,  which  is  ap¬ 
propriate  for  air.  In  the  first  experiment,  only  a  plasma  arc  was  accelerated, 
whereas  in  the  second  both  the  arc  and  a  2.5  g  projectile  were  employed.  The 

f* 

°K.A.  Jamison  and  H.S.  Burden  ( private  communication). 
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3.0 


Figure  2 .  The  parameters  and  z  ,  used  to  determine  the  characteristic  time 
and  the  time  of  maximum  velocity,  as  a  function  of  z^.[See  Eq .  (6).] 
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mass  of  the  arc  was  not  known  but  we  assumed  a  value  of  0.1  g,  as  has  been 

7  8 

done  in  previous  work.  *  Corresponding  values  of  zQ  for  the  two  experiments 
were  8.2  and  0.31,  respectively. 

Results  for  the  two  cases  are  shown  in  Fig.  3.  The  time  scale  is  in  mil¬ 
liseconds,  while  the  velocity  scale  is  in  km/s  for  Case  1  (arc  only,  no  pro¬ 
jectile)  and  hundreds  of  m/s  for  Case  2  (arc  and  projectile).  In  both  instances 
results  are  plotted  only  for  the  approximate  time  during  which  the  projectile 
or  armature  was  in  the  lm-long  gun  tube. 

As  is  evident  from  the  figure,  the  air  has  little  effect  upon  the  accel¬ 
eration  of  the  finite-sized  projectile  since  v^  differs  only  slightly  from  v 

during  the  course  of  acceleration.  In  fact,  the  characteristic  time  for  this 

case,  calculated  in  the  manner  described  previously,  is  1.25  ms  and  the  two 

velocities  do  not  differ  by  more  than  5%  for  times  less  than  t  .  It  should  be 

c 

noted  that  the  velocity  reaches  its  maximum  value  (^  631  m/s)  around  1.45  ms 
and  is  limited  primarily  by  the  time  decay  of  the  current.  The  atmosphere  has 
little  effect. 

For  the  case  in  which  only  the  arc  is  accelerated,  we  see  that  the  air 
exerts  a  very  significant  effect  upon  the  acceleration  process.  The  charac¬ 
teristic  time  for  this  case  is  only  0.04  ms  and  the  difference  between  v  and 
vq  grows  rapidly  thereafter.  The  maximum  velocity  reached  is  about  3.4  km/s 

and  is  attained  at  t  =  0.17  ms.  In  this  case  the  velocity  is  limited  primar¬ 
ily  by  the  effect  of  the  air  since  the  current  decays  little  during  the  time 
of  acceleration. 

Results  of  the  calculation  are  in  good  qualitative  agreement  with  the  ex¬ 
periments  in  question.  For  the  2.5  g  projectile,  the  velocity  increased 
steadily  over  the  first  37  cm  of  the  gun  tube  (the  location  of  the  last  detec¬ 
tor)  .  When  only  the  arc  was  accelerated,  however,  the  velocity  was  found  to 
rapidly  approach  about  3  km/s  and  remain  nearly  constant  thereafter,  behavior 
consistent  with  the  results  of  Fig.  3.  Specifically,  a  distance  vs.  time  curve 
plotted  at  the  location  of  the  four  detectors  (at  7,  17,  27,  and  37  cm  from 
the  breech  end  of  the  gun  tube)  produced  nearly  a  straight  line.  It  might  ap¬ 
pear  surprising  that  the  maximum  velocity  of  the  arc  is  predicted  reasonably 
accurately  since  the  arc  mass  was  only  crudely  estimated.  For  constant  cur¬ 
rents,  however,  which  nearly  apply  to  this  case,  v  becomes  independent  of 
m  and  is  given  by 


To  summarize,  the  nearly  constant  value  of  the  velocity  observed  in  the 
experiment  when  only  the  arc  was  accelerated  suggested  to  us  that  atmospheric 
drag  might  be  significant  in  some  rail -gun  experiments,  Results  of  the  calcu¬ 
lation  have  suggested  that  such  an  effect  is  indeed  important  at  high  velocities 

?I.E.  MoNab ,  J-  APPl-  Phys.  51,  2549  (1980). 

&J.D.  Powell  and  J.H.  Batteh,  J.  Appl.  Phys.  52,  2717  (1981). 
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t  (ms) 


Figure  3,  Arc  and  projectile  velocity,  with  and  without  atmospheric  effects, 

for  two  different  cases.  For  m  =  O.lg  (no  projectile),  the  vertical 
scale  is  in  km/s;  for  m  =  2t6g  (arc  and  projectile  present),  the 
vertical  scale  is  in  hundreds  of  m/s. 
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and  is  apparently  capable  of  explaining  these  experimental  results.  Precise 
conditions  under  which  air  exerts  an  appreciable  decelerating  force  on  the  pro¬ 
jectile  can  be  determined  for  any  experiment  by  calculating  zQ  for  the  case  at 

hand  and  using  Fig.  2  to  obtain  the  corresponding  characteristic  time. 
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ATTN:  Mr.  H.B.  Odom,  Code  F-12 
Dr,  M.F.  Rose,  Code  F-04 
Mr.  P.T.  Adams  and 

Mr.  D.L.  Brunson,  Code  G-35 
Dahlgren,  VA  22448 

2  Commander 

Naval  Research  Laboratory 
ATTN:  Dick  Ford,  Code  4774 
Dr.  I.  Vitkovitskv 
Washington,  DC  20375 

1  HQ  AFSC/XRB  (CPT  Dennis  Kirlin)/SDOA 
Andrews  AFB,  MD  20334 

2  AFATL/DLDB  (Lanny  Burdge,  Bill  Lucas) 
Eglin  AFB,  FL  32542 

1  AFATL  (Richard  Walley) 

Eglin  AFB,  FL  32542 

1  AFWL  (Dr.  William  L.  Baker) 

Kirtland  AFB,  NM  87117 

1  AFWL/NTYP  (John  Generosa) 

Kirtland  AFB,  NM  87117 

1  AFWL/SUL 

Kirtland  AFB,  NM  87117 

1  AFAPL  (Dr.  Charles  E.  Oberly) 

Wright -Patterson  AFB,  OH  45433 

1  AFAWL/P00S-2  (CPT  Jerry  Clark) 
Wright-Patterson  AFB,  OH  45433 

1  Director 

Brookhaven  National  Laboratory 
ATTN:  Dr.  James  R.  Powell,  Bldg  129 
Upton,  NY  11973 

2  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  R.S.  Hawke 

L-156,  Box  808 
Livermore,  CA  94550 
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DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


2  Director 

Los  Alamos  Scientific  Laboratory 
ATTN:  Dr.  Clarence  M.  Fowler,  MS970 
Dr.  Denis  R.  Peterson,  MS985 
Los  Alamos,  NM  87544 

3  NASA-LeRC-MS  501-7 
ATTN:  Bill  Kerslake 

Frank  Terdan 
Mike  Brasher 
2100  Brook  Park  RD 
Cleveland,  OH  44135 

1  Boeing  Aerospace  Company 
ATTN:  J.E.  Shrader 
PO  Box  3999 
Seattle,  WA  98124 

1  General  Dynamics 

ATTN:  Dr.  Jaime  Cuadros 
PO  Box  2507 
Pomona,  CA  91766 

3  GT  Devices 

ATTN:  Dr.  Derik  Tidman 

Dr.  Shyke  Goldstein 
Dr.  Neils  Winsor 
5705-A  General  Washington  Drive 
Alexandria,  VA  22312 

1  Physics  International 
ATTN:  Dr.  A.L.  Brooks 
2700  Merced  Street 
San  Leandro,  CA  94577 

2  R$D  Associates 
ATTN:  Mr.  Ronald  Cunningham 

Dr.  Peter  Turchi 
PO  Box  9695 

Marina  del  Rey,  CA  90291 

1  Science  Applications,  Inc. 

ATTN:  Dr.  Jad  H.  Batteh 
1503  Johnson  Ferry  RD,  Suite  100 
Marietta,  GA  30062 


Science  Applications,  Inc. 

Corporate  Headquarters 
ATTN:  Dr.  Frank  Chilton 
1250  Prospect  Plaza 

La  Jolla,  CA  92037 

Unidynamics/Phoenix,  Inc. 

ATTN:  W.R.  Richardson 
PO  Box  2990 
Phoenix,  AZ  85062 

2  Vought  Corporation 

ATTN:  William  B.  Freeman 
Charles  Haight 
PO  Box  225907 
Dallas,  TX  75265 

2  Westinghouse  Research  8  Development 
Laboratory 

ATTN :  Dr .  I an  R .  McNab 
Dr.  Y.  Thio 
1310  Beulah  RD 
Pittsburgh,  PA  15253 

2  Massachusetts  Institute  of  Technology 
Francis  Bitter  National  Magnet  Lab 
ATTN:  Dr.  Henry  H.  Kolm 
Dr.  Peter  Mongeau 
77  Massachusetts  Avenue 
Cambridge,  MA  02139 

1  University  of  Texas 

Center  of  Electromechanics 
ATTN:  Mr.  William  F.  Weldon 
167  Taylor  Hall 
Austin,  TX  78712 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 

ATTN:  DRXSY-D 

DRXSY-MP ,  H.  Cohen 

Cdr,  USATECOM 

ATTN:  DRSTE-TO-F 

Cdr,  CROC,  AMCCOM 

ATTN:  DRSMC-CLB-PA 
DRSMC-CLN 
DRSMC-CLJ-L 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


